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Using Monte Carlo simulations, we systematically investigate the non-equilibrium dynamics of 
the chiral degree of freedom in the two-dimensional fully frustrated XY model. The critical initial 
increase of the staggered chiral magnetization is observed. By means of the short-time dynamics 
approach, we estimate the second order phase transition temperature T c and all the dynamic and 
static critical exponents 0, z, (3 and v. 
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For critical phenomena, it is traditionally believed that 
universal behaviour exists only in equilibrium or in the 
long-time regime of the dynamical evolution. The univer- 
sal scaling behaviour is described by a number of critical 
exponents. Due to critical slowing down, numerical mea- 
surements of the critical exponents are very difficult. 

Recently, much progress has been achieved in dynamic 
critical phenomena. It was discovered that, already after 
a microscopic time scale t m i c , universal scaling behaviour 
emerges in the macroscopic short-time regime of the dy- 
namic process |IJ-^|. To illustrate this, we consider the 
following dynamic relaxation process: a magnetic system 
initially in a high-temperature state with a small initial 
magnetization too is quenched to the critical temperature 
T c without an external magnetic field and then released 
to a dynamic evolution with model A dynamics [Q. At 
the onset of the evolution, the magnetization is subjected 
to the scaling form |,|,H§ 



M(t, t, toq) ~ m t e F{t 1/uz T). 



(1) 



The exponent is a new independent exponent, t ~ 
(T — T c )/T c is the reduced temperature, (3 and v are the 
static critical exponents, and z is the dynamic exponent. 
At the critical temperature, r = 0, the magnetization 
undergoes a critical initial increase M(t) ~ t 6 '. 

In fact, the short-time dynamic scaling is very gen- 
eral. Another important example is the dynamic relax- 
ation of a magn etic system starting from an ordered state 
(too = 1) [fu^JIllll . The scaling form for this dynamic 
process is given by 

M^ k \t,T,L) = b~ k Pl v (b~ z t, b x l v T, b^L), (2) 

where, for later convenience, a system with finite size L 
has been considered. This finite size scaling has the same 
form as in the long-time regime but is assumed valid in 
the macroscopic short-time regime. 



One of the most prominent properties of the scaling 
forms (0) and (^|) is that the exponents (3, v and z take 
the same values as in equilibrium or in the long-time 
regime of the dynamic evolution. It has been suggested 
that it is possible to determine not only the dynamic but 
also all the static exponents as well as the critical temper- 
ature already in the short-time regime |]|| (see also Refs. 
jll] [HID). The method may be an alternative way for 
overcoming critical slowing down since the measurement 
does not enter the long-time regime. It is important to 
systematically verify this application in general and com- 
plex models. 

The two-dimensional fully frustrated XY (FFXY) 
model has been the topic of many recent studies [p^-p2[ . 
Critical properties of this model are rather unconven- 
tional. On the square lattice, the model has two kinds 
of phase transitions, i.e. the Kosterlitz-Thouless phase 
transition (XY-like) and the second order phase tran- 
sition (Ising-likc). Numerical simulations of the FFXY 
model suffer severely from critical slowing down. Due 
to the frustration, the standard cluster algorithm does 
not apply to the FFXY model. Most of the recent work 
on FFXY models supports that these two phase tran- 
sitions take place at two different temperatures, how- 
ever, their critical properties are still not very clear. For 
example, for the second order phase transition the es- 
timated values of the exponents (3 and v differ in the 
literature |2^,|2l],[l8|,[l7]] and the critical dynamics has 
not been investigated. It is still a matter of contro- 
versy whether the chiral degree of freedom of the FFXY 
model is in the same universality class as the Ising model 
p3|,p2|pll 



In this letter we present results of systematic Monte 
Carlo simulations for the short-time dynamic behaviour 
of the second order phase transition in the two- 
dimensional FFXY model. We confirm the short-time 
dynamic scaling. For the first time, we confidently de- 
termine the dynamic exponents and z. Based on the 
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short-time dynamic scaling, the static exponents and v 
as well as the critical temperature T c are also extracted 
from the numerical data. 

The Hamiltonian of the FFXY model on a square lat- 
tice can be written as 



H = —K 



<ij> 



cos 



(3) 



In our notation the factor 1/kT has been absorbed in 
the coupling K, Qi is the angle of the spin (a unit vector) 
located on site i, A^ determine the frustration and the 
sum is over the nearest neighbours. A simple realization 
of the FFXY model is by taking Aij = tt on half of the 
vertical links and = on other links. This is shown 
in Fig. ^ where the links with Ay = tt are marked by 
dotted lines. The spin configuration of one ground state 
of the model is also shown in Fig. |l|. The order parame- 
ter for the second order phase transition is the staggered 
chiral magnetization defined as (l9| 



sm(0i - 9j + A i:j ) 

<ij>£Pr 



(4) 



where (r x , r y ) is the coordinate of the plaquette P r . 



on the lattice, then we randomly choose a number of 
plaquettes and orient their spins according to the config- 
uration of the ground state shown in Fig. [I] till the initial 
magnetization mo is achieved. After the initial configu- 
ration is generated, the system is released to the dynamic 
evolution with the Metropolis algorithm at temperatures 
around T c . We have performed our simulation on lattices 
of size L = 128 and 256, and updated the system up to 
1000 Monte Carlo steps. The average is taken over in- 
dependent initial configurations with 40 000 samples for 
L = 128 and 10 000 samples for L = 256. Errors are 
estimated by dividing the samples into three groups. 
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FIG. 1. One of the ground states of the FFXY model. 



Dotted lines denote the links with Aij — tt. 



At first, we investigate the short-time critical be- 
haviour of Mj in the dynamic process starting from an 
initial state with a very high temperature and a small 
magnetization mo . We prepare an initial configuration in 
the following way: first we randomly generate all spins 



FIG. 2. The time evolution of the chiral magnetization 
Mi(t) when starting from a disordered initial state. The tem- 
peratures are T = 0.452, 0.454, and 0.456 from above. The 
dotted line represents the magnetization at T c = 0.4547. 

In order to locate the critical temperature T c , simu- 
lations have been carried out with three temperatures 
T = 0.452, 0.454 and 0.456. The initial magnetization 
is set to m = 0.06. In Fig. |[ the time evolution of the 
magnetization Mj (t) at different temperatures is plotted 
with solid lines in log- log scale for the lattice size L = 128. 
Data within the microscopic time scale t m i C ~ 100, which 
are dependent on microscopic details, are not included. 
Indeed we observe that the magnetization increases at 
the macroscopic early time. The magnetization at the 
temperature between T — 0.452 and 0.456 can be ob- 
tained by quadratic interpolation. From the scaling form 
(|l|) and as suggested in Ref. searching for a curve 
Mi(t) with the best power law behaviour can yield an 
estimate of the critical temperature. In Fig. ||, the dot- 
ted line represents such a curve and the corresponding 
critical temperature T c = 0.4547(8). From the slope of 
this curve, the exponent 8 = 0.200(3) is obtained. 

From an analysis of the data for L = 256 we have ob- 
served that the finite size effect for L = 128 is already 
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sufficiently small to be neglected. In fact, the finite size 
effect can easily be controlled in the short-time dynam- 
ics, and this is an advantage of the short-time dynamic 
approach. On the other hand, the critical exponent 9 is 
defined in the limit too — > 0. Therefore, the finite Too 
effect should also be considered. However, within sta- 
tistical errors our data for too = 0.04 and mo = 0.06 
show no difference. Hence the finite too effect will also 
be ignored. 

In principle, other static and dynamic critical expo- 
nents can be obtained from d T In Mj, the second mo- 
ment and the fourth moment of the magnetization, the 
auto-correlation and other observables. For example, 



\/vz= 0.59(3) is obtained from d T In Mj(t, r)| T= o. How- 
ever, this dynamic process starting from a disordered 
state is not the best choice to obtain these exponents or 
the critical temperature T c . A dynamic process starting 
from an ordered state is preferable, since the fluctuation 
is weaker. 

For this purpose, simulations were also performed with 
temperatures T = 0.452, 0.454 and 0.456, starting from 
an ordered initial state. The lattice size chosen was 
L = 256 and the system was updated for 2000 Monte 
Carlo steps. The average was taken over 2000 samples. 
The ordered initial state was taken to be the ground state 
shown in Fig. |l|. 
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FIG. 3. (a) The time evolutions of the chiral magnetizations Mi(t) when starting from an ordered state. From above, the 
solid lines represent the Mi(t) at T=0.452, 0.454, and 0.456. The dotted line is the magnetization at the critical temperature 
T c = 0.4545. (b) Deviation of Mi(t) from the power law behaviour within the time interval [200, 1000]. 



The estimation of T c can now be performed again. 
From the scaling form (||) and for sufficiently large L 
we can easily deduce the scaling behaviour for the mag- 
netization (k = 1) 

M/(t, r) = t- & t vz G{t x l vz T). (5) 

At the critical temperature, t = 0, Mj undergoes a power 
law decay. When r ^ this power law is modified by the 
scaling function Git 1 l vz T'). As pointed out earlier, the 
temperature for which the magnetization has the best 
power law behaviour is the critical temperature T c . In 
Fig. | (a), the time evolution of M>(t) at T = 0.452, 
0.454 and 0.456 is plotted in log-log scale. Mj(t) at 
other temperatures in the interval [0.452,0.456] can be 
estimated by a quadratic interpolation. To avoid the ef- 
fect of t m i c , measurements are performed in the time in- 
terval [200,2000]. The deviation from the power law can 
be estimated in different ways, e.g. as described in Ref. 



In this paper, we choose to measure the deviation 
as the error by fitting M/(i) directly to a power law in 
the time interval [200,2000]. Furthermore we perform 
the fitting in log scale, i.e. less weight is given to the 
data in the longer time regime. In Fig. || (b), the devi- 
ation of Mi from the power law is plotted as a function 
of the temperature. The clear minimum confidently in- 
dicates the critical temperature T c . The resulting value 
T c = 0.4545(2) is consistent with T c = 0.4547(8) ob- 
tained from Fig. | and very close to those values ranging 
from T c — 0.451 to 0.454 reported in most of the recent 
references p^G^pl]|22] . Our statistical error, however, is 
smaller. The corresponding magnetization is also plotted 
in Fig. H (a) with a dotted line. The slope of this curve 
yields the critical exponent (3/vz = 0.0602(2). The qual- 
ity of this measurement is very good. With [3/ v given, 
on can obtain a rigorous z or vice versa pfj| , p6| , pd] [9P] . As 
compared to simulations with a disordered initial state, 
the measurements here carry much less fluctuations. 
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FIG. 4. (a) The derivative d T In Mi(t, t)| t =o plotted versus time in log-log scale, (b) The time-dependent Binder cumulant 



U. 



To extract the critical exponent 1/vz, differentiation 
of Eq. (§) leads to 



9 T lnM J (t,r)| T=0 =t 1/,/ ^ lnG(r')| T '=o. 



(6) 



Therefore, d T In Mj(t, t)\ t=0 should also present a power 
law behaviour in the beginning of the time evolution. In 
Fig. | (a), d T In M/(t, r) at T c = 0.4545 is plotted in log- 
log scale. The power law behaviour is clearly seen. The 
slope yields the critical exponent 1/vz = 0.57(1). 

The final step is to seek the dynamical critical expo- 
nent z. For this, we introduce a time-dependent Binder 



cumulant U(t, L) = Mj 2 ^ /Mf — 1. Due to the short spa- 
tial correlation length in the short-time regime, a simple 
finite size scaling analysis shows at the critical tempera- 
ture 



U{t,L)~t d ' z 



(7) 



In Fig. [I] (b), the curve for U(t, L) in log-log scale shows 
a nice power law behaviour. The slope gives the critical 
exponent d/z = 0.92(2). Table [j] summarizes the results. 
The errors in the measurements from mo = 1 are clearly 
smaller than those from too = 0.06. 





T c 


e 


P/vz 


l/vz 


d/z 


7710 = 1 


0.4545(2) 




0.0602(2) 


0.57(1) 


0.92(2) 


mo = 0.06 


0.4547(8) 


0.200(3) 




0.59(3) 





TABLE I. The critical temperature and the critical exponents measured by simulations starting from an ordered as well as 
a disordered state. 





This work 


Ref. |2| 

(1996) 


Ref. y 
(1995) 


Ref. y 
(1994) 


Ref. y 
(1993) 


Ising 




0.4545(2) 


0.451(1) 


0.452(1) 


0.454(2) 


0.454(3) 




V 


0.81(2) 


0.898(3) 


1 


0.813(5) 


0.80(5) 


1 


2/3/iz 


0.261(5) 






0.22(2) 


0.38(2) 


0.25 


z 


2.17(4) 










2.165(10) 


e 


0.202(3) 










0.191(3) 



provide an estimate of the error on v = 1. For the Ising model, exponents v and 1/5 /v are exact values and 9 is taken from 
Refs. §,0. The exponent z in the literature ranges from 2.155 to 2.172 J6|p6|,[|,p7[ ■ Here an 'average' value is given. 
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The interesting and important property of the Binder 
cumulant is that one can estimate independently the dy- 
namic exponent z. With z in hand, we calculate the crit- 
ical exponents 2(3 /v and v from (3/vz and 1/vz. Table || 
lists all critical exponents along with the results reported 
in the recent literature. Now 9 — 0.202(3) is measured at 
T c = 0.4545, which shows a small difference from that at 
T = 0.4547. Our short-time dynamic measurements sup- 
port those from Ref . |Q and provide extra new results for 
the dynamic exponents z and 0. The exponent v of the 
FFXY model is different from that of the Ising model by 
nearly 20 per cent. This indicates that the chiral degree 
of freedom of the FFXY model is in a new universality 
class. Other exponents of the FFXY model do not differ 
much from those of the Ising model. From our numerical 
data, we observe that the values of the critical exponents 
are rather sensitive to the assumed or measured critical 
temperature T c . This should be one of the main reasons 
for the different values of the exponents reported in the 
literature. Figure || and, in particular, Fig. ||, give us 
confidence in our measurements of the critical tempera- 
ture T c . Furthermore, our exponent v is extracted from 
the data in the close neighbourhood of T c , in contrast 
to many simulations in equilibrium. The finite size ef- 
fect is also well under control in the short-time dynamic 
approach. 

In conclusions, using Monte Carlo methods, we have 
systematically investigated and confirmed the universal 
short-time dynamic behaviour of the second order phase 
transition in the two-dimensional fully frustrated XY 
model. Based on the short-time dynamic scaling form, all 
static and dynamical critical exponents are determined. 
The dynamic exponents 9 and z are obtained for the first 
time and the measurement of the exponent (3/vz and the 
critical temperature T c is very precise. The estimated 
value v = 0.81(2) is clearly different from v = 1 for the 
Ising model. Our investigation of the chiral degree of free- 
dom of the FFXY model is to date the most systematic. 
We are convinced that the short-time dynamic approach 
is not only conceptually interesting but also practically 
efficient. In the simulations we do not encounter diffi- 
culties associated with large correlation times since our 
measurements are carried out in the short-time regime 
and we do not have the problem of generating indepen- 
dent configurations. 

A possible extension of the present work is the 
Kosterlitz-Thouless phase transition. In fact, some ex- 
ponents have been obtained |2§| , |2"9f| . However, owing to 
the absence of symmetry breaking, a clear signal such as 
in Fig. H does not exist for the Kosterlitz-Thouless tran- 
sition temperature Tkt- The determination of Tkt & n d 
of the exponent v is very difficult and requires extensive 
simulations and a careful analysis. 
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